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Abstract 

The Bethe-Sommerfeld conjecture states that the spectrum of the stationary Schrodinger oper- 
ator with a periodic potential in dimensions higher than 1 has only finitely many gaps. After work 
done by many authors, it has been proven by now in full generality. Another case of a significant 
interest, due to its importance for the photonic crystal theory, is of a periodic Maxwell operator, 
where apparently no results of such kind are known. We establish here that in the case of a 2D 
photonic crystal, i.e. of the medium periodic in two variables and homogeneous in the third one, 
if the dielectric function is separable, the number of spectral gaps of the corresponding Maxwell 
operator is indeed finite. It is also shown that, as one would expect, when the medium is near to 
being homogeneous, there are no spectral gaps at all. 

1 Introduction 

The Bethe-Sommerfeld conjecture HI states that the spectrum of the stationary Schrodinger operator 

- A + V{x) (1) 

with a periodic potential V{x) in M", when n>2, has only finitely many gaps. Starting with [[T6l l2l 
and up to [|T3l . after work done by many authors, it has been proven by now in full generality (see 
[fTTl WSi ) for the history and detailed references). In presence of a periodic magnetic potential, the 
situation becomes much more complex. The corresponding result was proven in 2D case in [1 1] and 
jH. The proofs in the latter papers are very technical. In particular, [jTll used microlocal analysis 
tools of BH. Very recently L. Pamovski and A. Sobolev [[T4| have settled a much more general case, 
which allows in particular inclusion of magnetic terms. Another case of a significant interest is of the 
Maxwell operator in a periodic medium, where apparently no results of such kind are known. The 
importance of this problem stems from the photonic crystal theory (e.g., L2l|4l|5l[Tpj), where existence 
of spectral gaps is a major issue. We establish here that in the case of a 2D photonic crystal, i.e. of 
the medium periodic in two variables and homogeneous in the third one, if the dielectric function is 
separable, the number of spectral gaps of the corresponding Maxwell operator is indeed finite. It is 
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also shown that, as one would expect, when the medium is near to being homogeneous, there are no 

spectral gaps at all. 

Let us start with describing the relevant mathematical model. The standard form of the material 
Maxwell equations is 

'V - D = 47rp, 
V-B = 0, 

VxH=fJ + if . 

Here E and H are electric and magnetic fields, D is the electric displacement, B - the magnetic 
induction, p - the charge density, c - the speed of light, and J - the free current density. The fields 
E, H, B, D, and J are vector-valued functions from (or a subset of R^) into M^. We will assume 
absence of free charges and currents, that is p = and J = 0. 

We are interested in the EM wave propagation in an isotropic dielectric photonic crystal. In this 
case, the Maxwell equations should be supplemented by the constitutive (or material) equations 

D = £E, B = /xH. (3) 

Here e and ji are scalar time-independent functions called electric permittivity and magnetic perme- 
ability, correspondingly. In most photonic crystals considerations it is assumed that the material is 
nonmagnetic, that is /i = 1. We will also assume that the medium is periodic, that is e and ji are 
periodic with respect to a lattice F in M^. In what follows, we will assume F to coincide with the 
integer lattice 7?. 

Under the above assumptions, the Maxwell system reduces to the form 



VxE=-i^, V-H = 0, 
VxH = i£(a;)f, V • sE = 



^ (4) 



For mono-chromatic waves of frequency e M, one has E(a;,t) = e*'^*E(a;), H(x,t) = e*'^*H(a;), 
and thus one arrives to the spectral problem 



-'-V \ / E \ w / E 



^v^ y V H y c v H 



E 



(5) 



on the subspace S of smooth vector fields ^ jj J satisfying 

V • £E = 0, V • AtH = 0. (6) 

The operator 



-iV 



X 



is called the Maxwell operator. We consider M as the operator on the subspace S. 

We can extend the operator M to a self-adjoint operator acting on a Hilbert space %. The Hilbert 
space is a closed subspace of I-Lq — L2(R^, C^, sdx) ® L2(R^, C^,//dx). Namely, % consists 
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of those vector fields (E, H) G T-Lq for which V ■ e'E = 0, V ■ /iH = 0, where the divergence is 
understood in the distributional sense. Note that "H is the closure of 5 fl "Ho in the Hilbert space "Ho- 
The operator M is naturally extended to act on the set 

D = {(E,H) G H I V X E G L2(M^C^/idx), V x H G L2(M^ C^ edx)}. 

Here, as before, the differentiations are understood in the distributional sense. Now M is a self-adjoint 
operator on "H (see [5 J for more details). In what follows, we only need to know that co/c G M is in 
the spectrum of the Maxwell operator if the system ^ has a nonzero bounded solution (E, H) G S. 

One of the fields E or H could be eliminated and the problem can be re-written in terms of another. 
For instance, one can re-write dS]), © as the following second order spectral problem: 

'V X V X E = Xe{x)E 
V ■ eE = 0, 



where the spectral parameter A is equal to (7)^- Note that when e and fi are real- valued functions 
the spectrum of the Maxwell operator is symmetric with respect to the origin. Therefore A = (7)^ is 
in the spectrum of the generalized spectral problem ([8]) if and only if both uj/c and —u/c are in the 
spectrum of M (see ^ for the analogous conclusion). 

Our principal task is to show that under appropriate conditions on the periodic dielectric function, 
the spectrum of the problem dS]), and hence the spectrum of the operator M, has only finitely many 
gaps. While we expect this statement to hold in general, this text is devoted to proving it in a special 
case when 

e(a;i,X2,X3) = £:i(a;i) +£:2(a;2)- (9) 
We will show that the number of spectral gaps is finite, even if we restrict our consideration to the 
invariant subspace of the electric fields E that are normal to the plane (xi,X2) of periodicity and 
depend on (xi, X2) only. 

One should notice that, in spite of many similarities, there are some important differences be- 
tween the spectral problems for Schrodinger and Maxwell operators. This difference arises due to the 
multiplicative rather that additive appearance of the spectral parameter. This, in particular, applies to 
existence and location of gaps (see, e.g., [lOJ). It is easy to create gaps at the bottom of the spectrum 
of a periodic Schrodinger operator (for instance, creating a periodic array of well separated poten- 
tial wells, see [lOJ)- On the other hand, the spectrum of the problem ^ always starts at zero, thus 
preventing a similar gap opening approach. 

The paper is structured as follows: the main results (Theorems 12. 1[ 12.31 and 12.41) are stated in 
Section 2. It is also noticed there that Theorems 12.31 and 12.41 imply Theorem 12.11 In Section 3, the 
proof of Theorem 12.31 is reduced to an auxiliary Proposition 13.21 This Proposition, as well as other 
auxiliary statements, are proven in Sections 4 and 5. Section 6 contains the proof of Theorem 12.41 
The final Section 7 is devoted to final remarks and acknowledgments. 

2 Statement of the results 

Under the imposed assumption the Maxwell operator ^ admits an invariant subspace So C 
S of fields (E, H) = (Ei, E2, E3, Hi, H2, H^) that do not depend on X3. Furthermore, the space 
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is decomposed into the direct sum of two subspaces 5*1 © 5*2, where 5*1 consists of the fields 
{El, E2, 0, 0, 0, H) and 5*2 consists of the fields (0, 0, E, Hi, H2, 0). In physical terms, Si consists 
of the transverse electric (TE) polarized fields, while 5*2 consists of the transverse magnetic (TM) 
polarized fields. It is easy to observe that both Si and ^2 are invariant under the operator M. To show 
that the Maxwell operator has finitely many gaps or no gaps at all, it is enough to consider only TM 
polarized fields. In terms of the spectral problem ([8]), we assume that E = (0, 0, E{xi, X2)). Then the 
problem ^ reduces to the 2D scalar spectral problem 

-AE = Xe{x)E. (10) 
Thus the spectrum of problem ^ contains the spectrum of the operator — A considered as a self- 

E 

adjoint operator on the Hilbert space L2{M.'^, e{xi, X2) dxi dx2). 
Our main result is: 

Theorem 2.1. Let e{xi,X2,X'i) = ei{xi) + £2(3^2), where ei and 62 are C^-smooth, positive, 1- 
periodic functions on M. Then, 

1. The spectrum of the problem diOl) . and hence of the Maxwell operator M, contains a ray and 
thus has only finitely many gaps. 

2. If the functions Si and 62 are sufficiently close to constants uniformly on the whole real axis, 
then the spectrum o/470l) coincides with [0, 00) and has no gaps at all ( in this case, the spectrum 
of M coincides with the whole real axis). 

The well known Bloch theorem (see, e.g., its most general formulation in (9', Theorem 4.3.1]) 
provides a nice description of the spectrum of elliptic differential operators with periodic coefficients. 
In our case this theorem can be formulated as follows. 

Proposition 2.2. (e.g., /9) Theorem 4.3.1]) 

Let e G C^iM?) be a positive function periodic with respect to a lattice liL © I2L, i.e., 

e{xi + liui, X2 + 12^2) = e{xi,X2) 

for all Xi, X2 G M and ni,n2 G Z. Then the following statements are equivalent: 

L A number A > in the spectrum of the problem d70|) (in other, words, it is in the spectrum of 

the operator — Aj. 

e 

2. The differential equation 

-AE = \eE (11) 

has a bounded nonzero solution E. 

3. The equation 477]) has a nonzero Floquet-Bloch solution E that satisfies a cyclic (Floquet) 
condition 

E{xi + lini,X2 + /2^2) = E{xi, a;2)e^("-i+/3-2) 
for some a, /3 G M and all Xi,X2 G M, rii, n2 G N. 
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The last of the three statements above is the most convenient for us. We thus prove the following 
statement, which, according to Proposition l2.2l implies the first statement of Theorem 12.11 

Theorem 2.3. Let e{xi,X2) = Si{xi) + £2(3^2). where ei, 62 are C^-smooth, strictly positive, 1- 
periodic functions on M. Then there exists Aq > such that for any A > Aq the partial differential 
equation 

-AE = XeE 
has a bounded nonzero Floquet-Bloch solution 

E{xi,X2) = Ei{xi)E2{x2), 

where Ei, E2 are such that 

Ei{x, + 1) = e"'E,{xi),E2{x2 + 1) = 6^^^2(0:2), (12) 

with a, P eR. 

Furthermore, Aq, depends only on the number 

C:= . me.x mx)\,\e'Ax)\,\e':{x)l\{e.{x))-'\}. 

t=l,2,x€K 

The second statement of Theorem 12.11 follows if we establish the following result: 

Theorem 2.4. Let e G C(]R^) be a positive 1? -periodic function. Then, for any A > 0, there exists 
5 > such thatif\e{x) — 1| < 6, then for any < X < A, the partial differential equation 

-AE = XeE 

has a bounded nonzero solution in M^. 

The particular choice of the period (and thus lattice F) is not important for the proofs and can 
be made arbitrary by rescaling. For simplicity we will assume, as we have already agreed before, 
that F = Z^, and in particular "periodicity" of a function of one variable, unless specified otherwise, 
always means "1 -periodicity." 



3 Proof of Theorem 2.3 



We start with the standard separation of the variables and thus reduction to a one-dimensional prob- 
lem: 

Lemma 3.1. Let Si, £2 be continuous functions on M and A, c G M. Suppose that Ei is a solution to 
the differential equation 

E'l{x) + X{ei{x) + c)Ei{x) = 
and E2 is a solution to the differential equation 

E'^{x) + X{e2{x)-c)E2{x)=0. 
Then, the function E{xi,X2) = Ei{x)E2{x) is a solution to the partial differential equation 

-AE = XeE, 

where e{xi,X2) = ei{xi) + e2{x2). 
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The proof is straightforward. 

The proof of Theorem [2]3] will be extracted from the following auxiliary result: 

Proposition 3.2. Let e G C^(M) be a positive 1 -periodic function on M. Let do > be a constant 
such that 

max{|£:(x)|, \{e{x)y^\, \e'{x)\, \£"{x)\} < do. 

Then there exist positive constants Aq, c?i, d2 that depend only on do, such that the following property 
holds: 

If for some A > Aq the equation 

E" + XeE = Q 

does not have any bounded nonzero solutions, then the equation 

E" + X{e + c)E = Q 
has such a solution for any constant c satisfying 

-r < c < -!=■ 
A WA 

The proof of Proposition 13.21 will be provided in Section \5\ Now we are going to show how 
Theorem [23] can be derived from this proposition. 



3.1 Proof of Theorem 2.3 



Let us choose a constant do > such that 

2 d 

-r < £i < ^,Wi\ < do, W-l < do for i = 1, 2. 

Let also Aq, di, c?2 be constants provided by Proposition l3.2l for this particular do- We introduce a new 
constant 

L . . 2dr /2rfiY 
Ao = max <j Aq, cto«i, \ ~f~ j 

We will show now that any A greater than Aq is in the spectrum of (fTOl) . which will prove Theorem 

lO 

di di 
Let ci = 0, C2 = — , and C3 = — -. 

A A 

d 2d d 

We have \cj\ < — for j = 1, 2, 3. Since Aq > dodi and Aq > —z—, we obtain that \cj\ < — and 
^ Aq do 2 

also |c, | < — . It follows that 

do ^ 

— < £i + Cj < do 
do 
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for z = 1, 2 and j = 1, 2, 3. Besides, we clearly have 

\{ei + CjYl = \e'i\< do 

and 

\{ei + c,y'\ = \e'^\ < do. 

Hence the conclusion of Proposition 13.21 holds for either of the functions Si + Cj. 

Let J, A; G {1,2,3}, j ^ k. Then 



di 2di 
y<|c,-c,|< — 



(2,d \ ^ 
— i ) , we obtain that 
d2 J 

, 2di 2di d2 

\Ci — Ck\ < — — < — < —;=. 

Now Proposition 13 .21 implies that the equation 

E'l{x) + A(ei(x) + Cj)Ei{x) = (13) 

does not admit a bounded nonzero solution for at most one value of j = 1,2, 3. 
Similarly, the equation 

E'^ix) + X{e2{x) - Cj)E2{x) = (14) 

does not admit a bounded nonzero solution for at most one value of j = 1, 2, 3. Thus, for at least one 
j = 1, 2, 3 both equations (fT3l) and (fT4l) admit bounded nonzero solutions Ei and E2. Then, according 
to Lemma [3H, E{xi, X2) = Ei{xi)E2{x2) is a solution of the partial differential equation 

-/^E{xi,X2) = X{ei{xi) + e2{x2))E{xi,X2). 

Since E is clearly bounded and nonzero, Theorem |2.2| implies that A is in the spectrum of the operator 

A. This proves Theorem 12 .31 □ 

£1 + £2 

4 On the spectra of one-dimensional problems 

In order to prove Proposition 13 .21 we need to conduct an auxiliary study of the spectrum of the one- 
dimensional differential operator 

(15) 



u \dx^ 

where u and p are /-periodic functions and u > 0. First of all, the Bloch theorem in one-dimensional 
case implies the following result: 
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Theorem 4.1. (e.g., [9* Theorem 4.3.1]) Let u G C^(M) be a positive I -periodic function and p be 
piecewise continuous l-periodic function. Then the following are equivalent: 



1. A is in the spectrum of the operator — ( — p 



1 



2. The differential equation 

E"{x) + Xu{x)E{x) + p{x)E{x) = 
has a bounded nonzero solution E. 

3. For some a G M, the cyclic (Floquet) boundary value problem 

E"{x) + \u{x)E{x) + p{x)E{x) = 0, 
E{1) = E(0)e^", E'{1) = E'(0)e^" ^ ^ 

has a nonzero solution. 

Since the problem is now fomiulated on a finite interval, the spectrum of (fT6l) is no longer contin- 
uous. The following lemma is standard: 



Lemma 4.2. 1. The spectrum of problem Um is discrete and consists of a nondecreasing se- 
quence 

Ai < A2 < A3 < . . . 

such that Xn — )■ 00. 

2. The eigenvalues satisfy the variational principle 

. . r f'-f- (Pf) ■ f 

An = mf sup —— — , 

dim V = n / ^ 

where /i ■ /2 = fi{x)f2{x) dx is the E^-scalar product, is a space of functions f G 
H^[0, 1] such that f{l) = e*"/(0), and V is a vector subspace ofH^. 

3. Each eigenvalue Xn depends continuously on the parameter a G M, n = 1, 2, . . .. 

Indeed, due to ellipticity and compactness of the interval, the analytic Fredholm theorem (e.g., 
|l9l Theorem 1.6.16]) implies that the spectrum either coincides with the whole complex plane, or is 
discrete. Since obviously large negative values of A are not in the spectrum, the first statement of the 
lemma follows. The second and third statements are also straightforward. 

In view of Lemma 14^21 the range of A^ as a function of a is a closed interval J„, called the nth 
band of the spectrum, and the entire spectrum of the operator (fT5l) in (M) is the union of these bands 
for 77, = 1, 2, .... The neighboring bands are either adjacent, or else they are separated by a gap. 
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Now let us consider problem (fT6l) for a fixed value of a. The eigenvalues Ai, A2, ... (and the 
corresponding eigenfunctions) can be explicitly computed in the special case when u is constant and 
p = 0. To estimate A„ in a general case, we will reduce the considerations to this special case. The 
first step here is to apply the Liouville transformation (see [72^), which allows one to reduce the 
problem (fT6l ) to a similar problem, but now with a constant function u and a different value of the 
period /. This is done in Lemmas 14.31 - 14. 51 The function p is altered as well, and this is why it has 
been included in the problem (fT6l) in the first place (in the applications we will have p = 0). Then the 
influence of the function p on the spectrum is estimated using the variational principle (see Lemma 

Let u > 0,u G C^[0, /]. Let us define a function 



ax) = / V^^r) dr. 
Jo 

Then ^ G C^[0, 1] maps the interval [0, 1] homeomorphically onto the interval [0, A], where A = ^(/). 
We denote by z its inverse function, which is defined on [0, A]. 
Let 

9=^-^-^. (17) 

The function 9 is bounded and its upper bound can be easily estimated in terms of function u. 
The proof of the following lemma follows by a straightforward calculation. 

Lemma 4.3. Let E be a function on [0, /] and we introduce a new function 

F{y) = -^E{z{y)). 

Then, 

1. For any A G M, the function E is a solution of the differential equation 

E"{x) + \u{x)E{x) = 0, X G [0, /] 
if and only if the function F is a solution to the differential equation 

F"{y) + \F{y) + 9{z{y))F{y) =0, yE [0, A]. 

2. If additionally u satisfies the periodicity conditions u{l) = u{0), u'{l) = u'{0), then the function 
F satisfies the Floquet conditions F{A) = F(0)e*", F'{A) = F'(0)e*° if and only if E satisfies 
the similar Floquet conditions E {I) = E(0)e^°, = E'(0)e^° 

Now we are going to compare the spectrum of the problem (fT6l) in the case p = with the 
spectrum of an explicitly solvable problem. 
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Lemma 4.4. Letu G C^[0, /] be positive on [0, /] and satisfy (with its first derivative) periodic bound- 
ary conditions. Let also Ai < A2 < A3 < . . . Z?^ the eigenvalues of the problem 



E"{x) + \u{x)E{x) = 0, 
E{1) = E(0)e*",E'(/) = E'(0)e'" 

Further /ef AJ < A2 < A3 < . . . be the eigenvalues of the problem 

F"{x) + XF{x) = 0, 
F{A) = F(0)e*",F'(A) = F'(0)e^' 

where 



(18) 



(19) 



A= / VM(r)rfr. (20) 
Then |A„ — A* | < sup \0\for all n > 1, where 9 is defined in 4771) . 



Proof Let, as before, ^(x) = \/u{t) dr, x E [0, /] and let z be its inverse function. 
According to Lemma |431 a function E is the solution to the boundary problem (fTSl) if and only if 
the function F{y) = ^u{z{y))E{z{y)) is the solution to the boundary problem 

F"iy) + XFiy) + piy)Fiy) = 0, [0,^(0], 

F{A) = F{0)e'", F{A) = F(0)e'", ^ ^ 

where p{y) = 6{z{y)). Hence we have a one-to-one correspondence between eigenfunctions of the 
problems (fTSl ) and (|2T]) which is linear and preserves eigenvalues. Therefore, the sequence Ai < A2 < 
A3 < ... is also the spectrum of the problem (|2T]) . 
By the variational principle (Lemma |43I) . we have 

A„ = mf sup — — , 

dim V = n / ^ Q 

/'■/' 

a; = mf sup -p-r, 

dim V = n / ^ Q 

where /i • /2 = fi{x)f2{x) dx, is a space of functions / G i^'^[0, A] such that /(A) = e*"/(0). 
< sup IpI for any nonzero function / G i/^, it follows that 



Since 



|A„ - A*| < sup IpI 



for all n > 1. Clearly, sup |p| = sup 16*1. This proves the lemma. □ 
Let Ai < A2 < A3 < . . . be the eigenvalues of the problem (fTSi) . We recall that A„ is actually a 
continuous function of the parameter a. Note that the estimate on A„ obtained in Lemma 1441 does not 
depend on a. This allows us to estimate the entire band J„, the range of the function A„(a). 
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Lemma 4.5. Let C = sup \ 9\, where 6 is the function defined by f lTTI) . and A be the number defined 
by d20l) . Then, assuming that 

Tx\n-lf 



one has 



Moreover, assuming that 



one has the inclusion 



7r^(n-l)^ 



-c 



C Jn. 



IT 



C Ji 



Proof Consider the eigenvalue problem 



F"{y) + XF{y) = 0, 
F{A) = F(0)e*",F'(A) = F'{0)e' 



(22) 



Its eigenfunctions and eigenvalues are 



A* 



a + 2iTk 
A 



fk{y) = exp 



iy{a + 2TTk) 
A 



k ez. 



(23) 



Let A^, Ag, ... denote the above eigenvalues arranged in ascending order. Then for a E [0, vr) we have 



« — 27r 

A 



a + 27rrz 



a — 2'nn 
A 



Besides, A„(27r — a) = A„(q;) for all n > 1 and all a. 
Let J* be the range of A* as a function of a, i.e. 



j: 



' A^ 



In particular, we have A* (ai) = 7r^(n — 1)^^"^ and A* (02) = Tr^n^A"^ for ai = 0, 0^2 = tt if n is 
odd and ai = tt, 0:2 = if n is even. 

By LemmalMl |A„(ai) - A;(ai)| < C and |A„(a2) - A;(a2)| < C. Thus the first statement of 
the lemma follows. To prove the second statement, it remains to notice that always belongs to Ji, 
since Ai(0) = 0. □ 
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5 Proof of Proposition 1X2 



Before proceeding to the proof of Proposition 13 .21 we need two more technical statements. 

Suppose that some A > lies in a spectral gap of the operator —u~^cf/dx'^. We would like to 
know whether it is possible to modify function u slightly, so that A becomes covered by a spectral 
band. Lemma |43] tells us that A is close to a number vr^n^A^^, where n is an integer and A is defined 
in (|20|) . The next lemma explains how we need to change A to move all such numbers away from A. 
The parameter b in the lemma should be understood as tt^A^^. 

Lemma 5.1. Let Bq, bi > 0, c> 0. Suppose that 

I A — bri^l < c 
for some b^ < b < bi and integer n>l. Then 

I A — 6m^| > c 

for any integer m>l, provided that 

<\b-b\< — ^ 



A 12VA 
and A > max 

Proof Assume that b satisfies conditions of the lemma. We show first that |A — bm?] > c. 

Indeed, since A > 2c, the inequality I A — bn"^] < c implies that - < - — - < 2. Then 

2 bw^ 

\X — bn \ > \b — b\n — \X — bn \ > — —n — c> — c > 2c — c = c. 

A A 

Now we are going to show that A — b{n + 1)^ < — c and A — b{n — 1)^ > c. This will complete the 
proof. 

We have 

A-6(n+l)2 = X-b{n + l)^ + {b-b){n+lf < c-b{2n + l) + -^{n + lf < c-2bn+^ 



2bn^f\ bl'^n^ fbX y/haVXn'^ [bX 2 r— /M 2 r— 

' <c-\/ — H — < c- W — + - a/ooA < c- \/ — — + -a/ooA < 



sVA" V2 3A - V23^"- V2 3 



''''^1-7=2 J- 



By a hypothesis of the lemma y/boX > 60c, hence 



'Mi;-^)<c-60c(f--L|<-e. 
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It remains to prove that A — b(n — 1)^ > c. Since A > 2c it is no loss to assume that n > 2. We 
have 



\-b{n- if = \ - h{n - if + {h + h){n - if > -c + h{2n - 1) 



(,3/2 



12VA 



'n-lY. 



Since n > 2, we have in — if < and 2n — I > n, then 



-c + b{2n- I] 



.3/2 
l2^A 



^3/2 

'n - if > -c + bn > -c + 

l2^A ~ 



bny/X ^/boby/Xn^ 



12A 



> 



b\ 1 




6 



&nA 1 



c+ \hr- 77 V ^oA > -c + W ^ - - V^oA = -c + V^oA ^ - - > 



-c + 60c 



2 6 
1 1 

7i~ 6 



1 



6 



> 30c. 



This finishes the proof of Lemma l5.1[ 

Let Eo be a continuous positive function on [0, 1]. For any c > —mfso, let 



□ 



A{c) = / \/eo(xY+~c dx. 



As we know from Lemma |431 the quantity tt'^A{c) ^ is closely related to the location of the spectral 
gaps of the operator —SQ^cP'/dx^. We need to know how it depends on c. 

Lemma 5.2. Suppose that di < Eq < d2, where di and ^2 are positive constants. Then for any 
c G [-f, f] we have 



TT 

OS. 



c\ < 



TT 



TT 



< 



TT 



and 



A2(c) ^2(0) 

27r2 



2dV 



271^ TT 

< 



< 



3^2 ~ A2(c) - di 
Proof The function A is well defined on [—3-5 ^] ^i^d smooth. We have 



A'{c) 



dc 



'1 

dc 



(A/£:o(a;) + c) (ix 



2Weo{x) + c 



It follows that {2^d^Tc)-^ < A'{c) < {2^/drTc)-\ Also, VdTTc < A{c) < ^fd^c. 
Let us introduce a new function B[c) = ir'^A^'^^c). Then 



B'(c) 



n'A'jc) 
2^3 (c) ■ 



Therefore, 



TT 



4(^2 + C)2 



< \B'ic)\ < 



4(di + c)2- 
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As a consequence, 



15(c)- 5(0)1 < 



4(di + r 



dr 



< 



vr^lcl 



4:di{di + c)' 



Similarly, 

15(c)- 5(0)1 > 

In the case |c| < c/i/2, we obtain 

15(c)- 5(0)1 < 

and 



Ad2{d2 + c)' 



TT^Id vr^lcl 



c 









|5(c) -5(0)1 > — 
' ^ ^ ^ ^' - 4^2(^2 + f) " 

Also, in this case we have 

^/dj2 < ^/di + c < A{c) < Vc?2 + c < ^3^2/2, 

hence 

TTT < B{c) < — , 
c5a2 «i 

which proves the lemma. □ 
Now we are ready to prove the central Proposition 13 .21 

Proof of Proposition 13.21 Recall that the function e satisfies d^^ < £ < do, \e'\ < do, \e"\ < do- Let 
us introduce constants 

27r^ 5 
ao = ai = 27r^rfo, 0=^4 + ^O' 

^ 24 ,2„ , al^^ ^ 360062 

cti = ^airfo®' "2 = „ ^ ,2 ' Ao = max<^2e, ,{2dod2) >. 

vr^ 67r^cto I o-o J 

Suppose that for some A > Aq the differential equation 

E" + X{e + c)E = (24) 

does not admit a bounded nonzero solution for c = 0. We shall show that this equation does admit 
such a solution for any constant c satisfying 

rfi . I I ^ d2 

A ^A 
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First we introduce the following functions defined at least for c > —(Iq 



1. 



A{c) = / v^£(r) + c dr, 



and 



9{x,c) 



e"(x) 



16{e{x) + cy 4(e(x) + c)2' 
We shall need some estimates regarding A and 9. 

Take any c such that di/X < \c\ < d2/\/X. Notice that |c| < l/{2do) since A > Aq > {2dod2 
Then Lemma [S!2] implies that 



TT 



IcI < 



A2(c) A2(0) 



2^2 



< 



and 



IT 



Moreover, since e + c > e — |c| > \e'\ < d^, and \e"\ < d^, we have 



1^1 < 



'12 



16(e + c) 



4(£ + c) 



< 



+ 



dn 



e. 



We assumed that the equation (|24l) does not admit a nonzero bounded solution for c = 0. By 
Theorem 14. 1[ A is not an eigenvalue of the operator According to Lemma |431 this implies 

that 



A 



A(0)2 

for some integer n> 1. By the above we have 



< sup |^(x,o)| < e 

a-e[o,i] 



TT 



TT 



> 



77" , , _ -K^ di 4ai6 



Qdl 



c > 



Qdl A 



and 



TT 



TT 



A2( 



^2(0) 



< 



^ "0| I / ^ '^0 ^2 



C < 



Since A > max 1 29, ^^^}, it follows from Lemma O that 



3/2 



2 VA 12yA' 



A 



A2( 



> e 



for any integer m > 1. 

Since sup \9{x, c)| < 6, Lemma |431 implies that A is an eigenvalue of the operator — or, 
equivalently, the equation (f24l) has a bounded nonzero solution. □ 
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6 Proof of Theorem 1X4 



Let £ be a positive continuous function on periodic with respect to the integer lattice. Suppose that 

A > belongs to the spectrum of the operator — A or, equivalently, the differential equation 

e 

-AE = XeE 

has a bounded nonzero solution in M^. According to the Bloch theorem (Theorem l2.2l) . we can choose 
the bounded solution that satisfies a Floquet condition 

E{xi + hni, X2 + l2n2) = E{xi, X2)e'(""^+^"^) (25) 

for some a, /3 G M and all xi,X2 G M, ni,n2 E N. Then the function is a solution of the following 
boundary value problem in the unit square with quasiperiodic (or Floquet) boundary conditions: 

-AE = XeE, 

E{1, X2) = e-E(0, X2), X2) = e*"||(0, X2), (26) 

^(xi, 1) = e^^E{x,, 0), 1) = e'^|f(a;i, 0). 



Conversely, any solution of the boundary value problem (1261) can be extended to a solution of the 
equation —AE = XeE in the entire plane that satisfies the Floquet condition (|25l) . 

The next standard statement collects the information about the spectrum of the problem (|26l ) that 
we will need to prove Theorem 12 .41 

Proposition 6.1. [75]/) 

1. The spectrum of the problem ^26\) is discrete. Its eigenvalues form an nondecreasing sequence 

< Ai < A2 < A3 < . . . , Xn^ 00. 

2. Each eigenvalue A„ is a continuous function ofa,P G M. 

3. Dependence of the eigenvalue A„ = A„(a, (3; e) on the function e is monotone. Namely, ife < e 
everywhere in the unit square, then A„(a, /3; e) < A„(a, /3; e). 

The next lemma provides, also a standard, statement on dependence on the dielectric function e. 

Lemma 6.2. For any d > and n there exists 5 > such that 

\Xn{a,(5]e) - Xn{a,l3]e)\ < d 
for all a, (3 provided that \e — e\ < 6. 
Proof. Clearly, 

Xnia,/3; ke) = -Xn{a, /3; e) 
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for any k > 0. Let us pick k > such that KA„(a, /3; e) < i!). Then 

A„(a, (3; (1 + > A„(a, /3; £:)>(!- K)A„(a, (3; e) > A„(a, /3; e) - {}. 

I + K 



Also 

A„ I rT. 

1 + K 

If the function e satisfies inequalities 

e 



A„ I a, (3;-^ ) < (1 + K)A„(a, /3; < A„(a, /3; e) + 



< £ < (1 + k)£ 



1 + K 

everywhere in the unit square, then, by Proposition 16. 11 



A„(a, (3; (1 + K)e) < A„(a, e) < A„ ( a, /3 ^ 



' 1 + K 



which implies that 

|A„(a,/3;e) - A„(a,/3;e)| < i9. 

Since 

inf £ = /i > 0, 

we have 

(1 + K)e — e = ne > Kfx, 
e-- = > ^ 



Thus, for any function e such that 

\s — e\ < < Ku, 

1 + K 

we have | A„(a;, (3; e) — A„(a;, (3;e)\ < i!). This proves the statement of the lemma. □ 
Let us fix the function e. For any n let J„ = /^(e) denote the range of A„(a;, /3; e) as a function of 
a and /9. It follows from Proposition 16. II that /„ is a closed interval. Note that this interval lies in the 

spectrum of the operator — A acting on the entire plane. Furthermore, the spectrum is exactly the 

6 

union of the intervals h{e), hi^), 

The following statement about the spectrum of the Laplace operator is well known (e.g., ifTTlD . 
and can be proven easily, so we skip its proof: 

Lemma 6.3. In the case e = 1, the intervals I\, I2, h,..- overlap. That is, for any n the intersection 
In n In+i has a nonempty interior. 
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Proof of Theorem |2.4i We can now address the proof of the remaining result, Theorem |2.4l which 
follows rather immediately from Lemmas 16.21 and [63l Indeed: 

The spectrum of the operator —e^^A is the union of the intervals Ii{e), hi^), ■ ■ ■ defined earlier in 
this section. Hence we need to show that for any A > this union covers the interval [0, A] provided 
that the function e is close enough to 1 uniformly. 

Let = [ai, 12(1) = [02, &2], According to Lemma [631 a^+i < 6„. Given A > 0, let 

us take such that > A. Let 

'&o = ^ min {bn - ctn+i)- 

2 l<n<N 

By definition, 2^9o is a lower bound on the length of the intersection /„ fl In+i for n = 1,2, . . . , N — 1. 
Further, let 1!) = mm{'dQ, A — b^}. 

We have a„, = A„(a„, 1) and 6„, = A„(a^, 1) for some «„, a'^, G M. By Lemma [Ol 
there exists 5 > such that 

|A„(a„,/3„;£) - a„| < {}, 

forn = 1, 2, . . . , A^ whenever the function e satisfies sup — 1 1 < 5. Since the points A„(a„, e) 
and A„(a;^, e) lie in the interval /„(£:), we obtain /.„(£:) D [a„ + 6„ — Moreover, D 
[0, 61 — ^9] as Ai(0, 0; £:) = (indeed, the constant function is an eigenfunction of the problem (|26l ) 
with periodic boundary conditions for A = 0). By the choice of 'd, the intervals and In+i{s) 
overlap for n = 1, 2, . . . , A^ — 1. Besides, the right end of the interval /Ar(e) lies to the right of 
the point A. Thus the intervals I2{s),...Jn{£) cover the interval [0, A] without any gaps. This 

finishes the proof of Theorem |2.4[ □ 

Final remarks and acknowledgments 

In this paper, we only considered the E-polarized modes for electromagnetic waves propagating along 
the periodicity plane of a 2D photonic crystal with a separable dielectric function. There clearly 
remain several issues to consider. The separability condition is a strong restriction, and thus one would 
want to avoid it. Besides, the case of fully 3D periodic photonic crystals has not been considered. 
The author plans to address these questions in the future work. 

The author is grateful to Peter Kuchment and Yaroslav Vorobets for helpful discussions. 
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